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ABSTRACT 

A  validation  is  carried  out  of  the  steady  state  thermal  stress  analysis  capability  of  the 
P AFEC  Finite  Element  program.  Two  problems  are  considered  for  which  closed  form 
solutions  can  be  obtained:  these  comprise  the  thermal  stresses  in  a  heated  circular  plate 
and  also  a  bonded  joint.  Also,  a  procedure  is  developed  for  calculation  of  residual 
stresses  in  bonded  joints.  Comparisons  with  closed  form  solutions  indicate  the  degree 
of  discretisation  necessary  to  achieve  accurate  results. 
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Development  and  Validation  of  a  Finite  Element 
Based  Method  to  Determine  Thermally  Induced 
Stresses  in  Bonded  Joints  of  Dissimilar  Materials 


Executive  Summary 

Acoirate  computational  thermal  stress  analysis  is  particularly  important  for 
assessment  of  Ae  viability  of  bonded  repairs  to  military  aircraft.  Thermal  stresses  arise 
due  to  the  different  coefficients  of  thermal  expansion  of  the  materials  involved  and 
may  be  significant  at  the  extremes  of  the  operating  temperatures  of  these  aircraft. 
Furthermore  it  is  known  that  the  residual  thermal  stresses  due  to  the  bonding  process 
may  be  significant.  Previously,  the  thermal  stress  capability  of  the  PAFEC  Finite 
Element  program  has  not  been  used  at  AMRL.  Also  no  PAFEC  steady  state 
temperature  test  problems  are  available.  As  a  result  it  has  been  necessary  to  carry  out 
teste  for  which  closed  form  solutions  can  be  obtained.  The  problems  considered  are  the 
thermal  stresses  in  a  heated  circular  plate  and  a  bonded  joint.  Furthermore  an 
algorithm  is  developed  for  the  calculation  of  residual  stress.  Results  obtained  show 
that  the  procedures  developed  here  give  an  accurate  assessment  of  the  thermally 
induced  stress  state  in  bonded  joints.  The  algorithm  developed  can  be  applied  to  more 
complex  geometries  to  enable  a  full  assessment  of  the  thermal  stress  aspects  of  bonded 
repairs. 
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Notation 


k  ply  ntimber  in  laminate 

n  total  number  of  plies  in  laminate 

r  radius 

s  stiffness  ratio 

t  thickness  of  laminate 

to  thickness  of  inner  region  of  circular  plate 

ti  thickness  of  plate 

t2  thickness  of  patch 

ta  thickness  of  adhesive 

Uo  displacement  of  inner  region  of  circular  plate 
ui  displacement  in  outer  region  in  circular  plate 
ui  displacement  in  ID  bonded  joint 
U2  displacement  in  patch  on  ciroilar  plate 
U2  displacement  of  reinforcement  in  ID  bonded  joint 

X  Cartesian  co-ordinate 

y  Cartesian  co-ordinate 

z  Cartesian  co-ordinate 

Cl  constant  of  integration 
Cz  constant  of  integration 
C3  constant  of  integration 
E  Yotmg's  modulus 

Eo  Young's  modulus  for  inner  region  of  circular  plate 

El  Yotmg's  modulus  in  circular  plate 

E2  Young's  modulus  in  patch  on  circular  plate 

Ell  Young's  modulus  ortitotropic  material  in  1  direction 

E22  Young's  modulus  orthotropic  material  in  2  direction 

G  shear  modulus 

Ga  shear  modulus  of  adhesive 

Ri  inner  radius  of  circular  plate 

Ro  outer  radius  of  circular  plate 

Ti  temperature  of  inner  region 

To  temperature  at  outer  radius 


Notation  continued 


a  coefficient  of  thermal  expansion 

ao  coefficient  of  thermal  expansion  of  inner  region  of  circular  plate 

tti  coefficient  of  thermal  expansion  of  circular  plate 

tti  coefficient  of  thermal  expansion  of  component  1  in  ID  bonded  joint 

0.2  coefficient  of  thermal  expansion  of  patch  on  circular  plate 

0.2  coefficient  of  thermal  expansion  of  reinforcement  in  ID  bonded  joint 

ttL  longitudinal  coefficient  of  thermal  expansion  of  laminate 

ot  transverse  coefficient  of  thermal  expansion  of  laminate 

Ueff  effective  coefficient  of  thermal  expansion  of  circular  plate 

V  Poisson's  ratio  in  isotropic  material 

vi2  Poisson's  ratio  for  contraction  in  2  direction  associated  with  extension  in  1 
direction 

V21  Poisson's  ratio  for  contraction  in  1  direction  associated  with  extension  in  2 
direction 
a  direct  stress 

Ca  appHed  stress  to  bonded  joint 

Go  radial  stress  in  inner  region  of  ciroilar  plate 

Gi  radial  stress  in  circular  plate 

Gi  direct  stress  in  ID  bonded  joint 

G2  radial  stress  in  patch  on  circular  plate 

G2  direct  stress  in  reinforcement  of  ID  bonded  joint 

Sr  radial  strain 

T  shear  stress  in  ID  bonded  joint 

Tmax  maximum  shear  stress  in  ID  bonded  joint 

0  angle  of  k*  ply  to  longitudinal  axes 

[B]  relates  strains  to  displacements 

[C]  conductivity  matrix 

[D]  elasticity  matrix 
[K]  stiffness  matrix 

{d}  displacement  vector 

{P}  load  vector 

{Pt}  thermal  load  vector 

{Q}  heat  flux  vector 

{T}  nodal  temperature  vector 

V2  Laplacian  operator 
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1.  Introduction 

In  this  report  comparisons  are  made  between  thie  F.E.  results  for  a  number  of  thermal 
stress  solutions  and  closed  form  solutions.  The  thermal  stress  options  for  the  PAFEC 
F.E.  program  have  not  been  previously  used  at  AMRL  and  a  validation  of  these 
solutions  is  required.  While  standard  test  data  files  are  available  for  transient  analysis, 
no  steady  state  only  test  problems  exist.  To  gain  confidence  in  the  results  it  is  necessary 
to  validate  the  solutions.  The  most  useful  solutions  are  those  for  which  differential 
equation  solutions  exist.  In  this  paper  the  heated  circular  plate  is  analysed  since  the 
closed  form  solution  is  available,  and  also  this  will  validate  a  two  dimensional 
analysis  The  solution  for  the  adhesive  shear  stress  in  a  bonded  joint  will  also  be 
considered.  Furthermore  the  F.E.  method  will  be  appHed  to  a  residual  stress  problem 
for  which  an  analytic  solution  is  available. 


2.  Finite  Element  Thennal  Stress  Analysis 


A  thermal  stress  analysis  is  usually  carried  out  with  the  intention  to  calculate 
thermally  induced  stresses,  strains  or  displacements.  Thermal  stresses  may  mise,  for 
example,  in  a  bonded  joint  consisting  of  materials  with  different  thermal  coefficients  of 
thermal  expansion.  Thermal  stresses  may  occur  in  a  heated  structure  which  is  rigidly 
constrained,  and  also  in  a  structure  with  temperature  gradients.  Thermal  stress 
analysis  is  usuaUy  carried  out  in  two  steps,  the  first  being  the  thermal  analysis  which 
will  calculate  temperatures  at  each  node,  and  die  second  to  calculate  the  corresponding 
stresses  and  displacements. 

The  simplest  thermal  analysis  is  that  of  steady  state  temperature  or  steady  state 
conduction.  In  this  case  the  equations  are  linear  and  are  derived  from  the  Laplace 
relationship  and  in  matrix  form  are: 

[Cir)  =  {Q]  -m 

where: 

[C]  is  the  conductivity  matrix 
{ T]  is  the  nodal  temperature  vector 
{Q]  is  the  heat  flux  vector. 

The  boundary  conditions  may  be  specified  as  temperatures  or  heat  flux  at  nodes.  A 
partitioning  of  equation  (1)  is  carried  out  into  unknown  and  known  temperatures.  If 
the  heat  flux  is  known  at  a  node,  then  the  temperature  is  treated  as  being  unknown  at 
that  node.  The  thermal  solution  carried  out  by  PAFEC  uses  'thermal'  elements  whose 
properties  are  dependent  on  the  geometry  and  tiiermal  conductivity  for  that  material. 
This  analysis  solves  for  the  temperature  vector  and  writes  it  to  a  file  for  use  in  the 
structural  analysis  which  follows. 
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The  Static  structural  solution  involves  the  matrix  equations: 


[/:]{4  =  {p}  +  {p,} 

where: 

[K] 

is  the  structure  stiffness  matrix 

W} 

is  the  displacement  vector 

{P] 

is  the  applied  load 

[Pt] 

is  the  thermal  load  vector  and  is  given  by: 

where: 

[5] 

relates  strains  to  displacements 

P] 

is  the  elasticity  matrix 

H 

are  the  coefficients  of  expansion 

[r] 

is  the  temperature  vector  from  the  thermal  solution 

dv 

integration  is  taken  over  the  voltime. 

The  thermal  load  vector  is  added  to  any  appHed  loads  that  exist,  and  the  usual 
displacement  solution  is  obtained.  This  corresponds  to  the  second  PAFEC  run  in  which 
the  previously  calculated  temperatures  are  read  back  in. 


3.  Heating  of  a  circular  plate 


Consider  now  a  circular  plate  shown  in  Fig.  1.  This  plate  is  uniformly  heated  within  a 
radius  r  =  R,  equal  to  a  temperature  Tj,  while  the  temperature  at  the  boundary  is  7^ . 
The  temperature  solution  given  in  ref(l)  which  satisfies  the  Laplacian  operator: 


is  given  by: 


v"r  =  o 


...(4) 


T=To  + 


(T,-To)ln(r/Ro) 

ln(R,/Ro) 

T=Ti 


for  Rj  <r<RQ  and 
for  r<Rj 


...(5) 


For  this  temperature  distribution,  Ihe  following  differential  equation,  ref(2),  is 
applicable  to  plane  stress: 
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where; 


dr  \r)  dr 


/I 

=  «(l+v)- 


r  is  the  radius 

u  is  the  radial  displacement 

V  is  Poisson's  ratio 

a  is  the  coefficient  of  thermal  expansion 

T  is  the  temperature. 


For  the  case  of  plane  stress,  and  fixed  edges  the  strain  can  be  defined  as; 

-  To) 


...(6) 


•(7) 


Where  a  ^  is  the  effective  coefficient  of  thermal  expansion  and  is  given  by; 


a(l+v) 


Ro(T,-To) 


RaT-(^)\rTdr-(^)irTdr 
r  0  Kn  0 


...(8) 


The  explicit  solution  given  by  ref(l)  for  the  case  of  fixed  edges  for  r<Rj  is: 


a(l  +  v)L^  l-Rj/Rl 

2  2ln(Rj  / Rq) 


...(9) 


As  expected,  the  expression  is  dependent  on  the  geometry  only  and  not  temperature. 

Firstly  a  heat  flow  analysis  is  carried  out  using  'thermal'  elements,  ref(3).  For  the 
example  considered  here  die  boxmdary  conditions  are  specified  as  Ti  =  lOCPC  at  r<R^ 
and  To  =  O^C  at  r  =  Ro.  Also  inner  radius  Ri  =  162  mm  and  outer  radius  Ro  =  5000  mm. 
The  F.E.  mesh  for  this  case  is  shown  in  Fig.2.  and  indicates  the  discretisation  necessary 
to  give  good  results.  Material  and  thermal  constants  are  given  in  Tables  1  and  2. 
Temperatures  computed  from  the  F.E.  analysis,  in  the  radial  direction,  are  presented  in 
Table  3  togedier  with  restdts  predicted  from  equation  (5).  In  this  case  agreement  has 
been  achieved  to  at  least  four  significant  figures  and  is  a  validation  of  this  part  of  the 
analysis.  Also  contour  plots  of  temperature  show  concentric  rings  about  r  =  0  and 
coirfirm  that  the  solution  is  radially  dependent  only. 

The  second  analysis  is  carried  out  with  the  geometry  above  but  with  equivalent 
structural  elements,  and  in  this  case  the  boundary  conditions  correspond  to  restraints 
around  file  outer  circumference  of  the  plate.  A  comparison  between  the  F.E.  and 
analytical  results  for  Ueff  are  shown  in  Table  4  for  a  range  of  outer  radii  (Ro).  The  F.E. 
values  of  Ueff  have  been  computed  from  displacements.  This  is  necessary  since  in 
PAFEC,  strains  are  calculated  from  stresses,  and  in  the  case  of  thermal  stresses  in 
which  the  structure  is  restrained,  these  strains  do  not  correspond  to  real  strains. 
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The  analytic  results  in  Table  4  are  calculated  from  equation  (9)  and  correspond  to  fixed 
edge  conditions  around  the  circumference.  The  comparison  between  F.E.  and  analytic 
results  shows  that  good  agreement  has  been  obtained  for  0  <r<R/.  Also,  contour  plots 
of  Von  Mises  strains  show  concentric  rings  about  r  =  0  and  cilso  indicate  that  the 
solution  is  radially  dependent  only. 

The  validation  obtained  here  is  particularly  important.  As  mentioned  in  ref(l)  die 
calculation  of  the  effective  coefficient  of  expansion  in  a  structure  is  important  in 
determining  the  residual  stresses  in  a  bonded  structure.  The  effective  coefficient  of 
expansion  of  the  plate  structure  is  very  much  lower  than  the  value  a  in  the  asymptotic 
limit  Rj  I  Rq  =  0, 


i.e. 


a(l+v) 


...(10) 


From  equation  (10)  the  asymptotic  value  for  the  plate  is  14.95x10^  and  the  maximum 
value  of  aeff  in  Table  4, 12.7xlOA  corresponding  to  the  maximum  radius,  is  stiU  less 
than  this  value. 

The  lower  value  of  the  coefficient  of  thermal  expansion,  in  this  example,  is  firstly  due 
to  the  circumferential  restraints  applied  to  the  edge  of  the  plate  and  secondly  is  the 
result  of  the  resistance  of  the  surrounding  colder  structure.  This  occurs,  for  example,  in 
the  bonding  of  a  high  strength  repair  to  the  wing  of  an  aircraft  in  which  a  high 
temperature,  over  a  local  area,  is  required  to  ciue  the  adhesive. 


4.  Thermal  stresses  in  a  bonded  joint 


4.1  1  Dimensional  equation 

In  a  bonded  joint  in  which  two  materials  are  bonded  together,  thermal  stresses  may 
develop  as  a  result  of  the  difference  in  thermal  coefficients  of  expansion.  The  following 
derivation  will  be  given  for  a  simple  double  overlap  joint  whose  geometry  is  given  in 
Fig.3  together  with  the  location  of  the  origin  of  the  x  axis.  It  is  assumed  that  aU  the  load 
transferred  by  the  lap  joint  is  by  adhesive  shear.  Also,  in  this  symmetrical  joint  it  is 
assumed  that  no  bending  takes  place. 

The  constitutive  equations  are: 


du,  <7 

— 5-  =  —  +  «,Ar 

dx  E,  ' 


dWj  O' 2  .rri 

— -  =  —^  +  a.,AT 
dx  E2 


...(11) 


where: 


M,  and  U2  are  the  displacements  in  the  components  of  the  joint 
tti  and  ai  are  the  thermal  coefficients  of  expansion 


4 


and  £3  Voting's  moduli  of  the  two  materials 

AT  is  the  change  of  temperature 

oi  and  02  are  the  stresses  in  the  components  of  the  joint. 
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The  equilibrium  equations  are: 
dx  f, 


^__L  =  0  ...(12) 

dx 


where: 


are  the  different  thicknesses  of  tiie  two  components  and 
T  is  the  shear  stress  in  the  adhesive. 


Compatibility  requires  that  the  shear  stress  is  given  by: 


where: 


is  the  thickness  of  the  adhesive 
G  is  the  shear  modulus  of  the  adhesive. 


...(13) 


Differentiating  equation  (13)  twice,  equation  (11)  once,  and  using  equation  (12)  leads  to 
a  simple  d.e.: 


dx^  _  £"2^2  -^1^1  _ 


...(14) 


It  can  be  shown  that  the  solution  of  this  equation  is  given  by: 


where: 


^a(^d\  ^2^2) 


Differentiating  equation  (13),  equation  (15)  and  using  equation  (11)  leads  to: 


h  ^1 


For  thermal  stresses  only  cr^  =  0  and: 
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r  =  ^^(ay-a2)ATe -(17) 

K 

The  direct  stress  for  component  1  is  determined  from  equation  (12)  and  equation  (17), 
and  is  given  by: 


GJ^ 

(j.  = - 


-^  +  (a,  -a2)Ar 

■^1 


...(18) 


At  x=0,  cr j  =  0 ,  and  rises  to  a  maximum  value  when  x  is  large  which  is  given  by: 


Gf 


o,  =-■ 


^  +  (a.-a,)Ar 


.(19) 


If  we  consider  a  single  lap  joint  in  which  bending  is  restrained  and  if  the  thickness  of 
the  skin  is  ,  then  =  tj  2. 


4.2  2D  F.E.  analysis  of  a  bonded  joint 

A  2D  F.E.  (plane  stress)  analysis  has  been  carried  out  for  a  bonded  joint  whose 
geometry  is  shown  in  Fig.4.  Only  2D  isotropic,  8  noded  quadrilateral  elements  are  used 
in  the  analysis,  and  aU  lie  in  the  xy  plane.  Fig.4  is  a  diagram  only,  but  does  indicate  the 
areas  in  wliich  the  mesh  is  refined.  Since  a  comparison  is  being  made  with  the  ID  d.e. 
which  was  developed  earlier,  the  mesh  contains  only  one  element  per  structural 
component  in  the  y  direction. 

In  the  thermal  analysis  the  temperature  has  been  set  at  a  uniform  temperature  of  lOQoC 
throughout  the  joint.  This  is  achieved  efficiently  by  specifying  the  temperature  at 
botmdaries  in  such  a  way  that  the  F.E.  program  completes  the  analysis  by 
interpolation.  The  relevant  Aermal  properties  are  coefficients  of  thermal  expansion, 
«,and  a^,  for  components  1  and  2  given  in  Table  5.  While  the  adhesive  also  has  a 
thermal  coefficient  of  expansion,  this  is  set  to  zero  to  enable  the  best  comparison  with 
the  ID  d.e.  Also,  it  has  been  foxmd  that  for  PAFEC,  non-zero  conductivities  must  be  set 
for  isotropic  elements.  In  steady  state  problems  the  actual  value  of  the  conductivity  has 
no  influence  on  the  results  provided  it  is  not  zero. 

In  the  structural  analysis  that  follows  the  corresponding  'structural'  elements  are  used 
and  the  material  properties  used  are  given  in  Table  5.  In  the  analysis  the  restraints 
used  are: 

1)  restraints  are  used  in  the  x  direction  for  all  nodes  located  on  component  1  with 
co-ordinates  of  x  =  0 
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2)  restraint  in  the  y  direction  for  all  nodes  in  component  1  with  the  co-ordinate  of 
y  =  0. 

The  restraints  in  the  y  direction  have  been  applied  to  prevent  bending  of  component  1, 
so  that  a  direct  comparison  can  be  made  with  the  ID  d.e..  Also  in  this  analysis  no  non- 
thermal  loads  have  been  applied. 


From  the  results  a  comparison  of  shear  stresses  can  be  made  in  two  ways.  Firstly,  the 
shear  stresses  are  caloilated  using  the  displacements  from  the  F.E.  analysis  together 
with  equation  (13).  Secondly,  the  shear  stress  is  taken  directly  from  the  F.E.  results 
corresponding  to  the  maximum  value  that  occurs  in  the  midplane  of  the  adhesive. 
These  results  are  shown  in  Table  6,  for  a  range  of  E2  /  B,  values,  together  with  the 
values  predicted  by  the  ID  d.e.  (It  should  be  noted  that  while  the  ratio  of  E2  /  E,  is 
being  varied  all  other  parameters  remain  the  same.)  As  shown  in  Table  6,  the  midplane 
solution  gives  better  agreement  with  the  ID  d.e.  than  the  solution  obtained  from 
displacements. 

Overall  the  difference  in  results  between  the  2D  F.E.  and  ID  d.e.  do  not  necessarily 
indicate  the  existence  of  an  error  in  either  solution.  Factors  that  may  influence  results 
are  firstly  the  Poisson's  ratio  effect  which  is  not  considered  in  the  ID  d.e.  Secondly,  in 
the  F.E.  analysis  shear  deformation  occurs  in  aU  components  of  the  joint,  while  in  the 
ID  d.e.  it  only  occurs  in  the  adhesive.  Furthermore,  the  F.E.  analysis  of  bonded  repairs 
in  general  shows  a  considerable  variation  in  shear  stress  from  the  patch/ adhesive 
interface  to  the  plate/ adhesive  interface. 

To  achieve  good  results  in  the  F.E.  analysis  a  very  fine  mesh  of  0.05mm  increments  has 
been  required  to  pick  up  dre  rise  from  x  =  0  to  a  maximum  at  up  to  0.25  mm  away  from 
the  end  of  the  joint.  In  the  case  of  the  d.e.  the  formulation  is  such  that  the  maximum 
value  occurs  at  the  end  of  the  joint.  Although  results  are  not  given  here,  the  decay  of 
the  F.E.  and  d.e.  differ  slightly. 

In  the  case  just  considered  the  F.E.  analysis  was  confined  to  the  analysis  of  isotropic 
materials.  Now  consider  the  case  in  which  component  2  has  orthotropic  properties. 
Composite  materials  are  orthotropic  and  bonded  joints  are  often  comprised  of  such 
materials.  In  the  repair  of  cracked  metallic  structures,  materials  such  as  boron/ epoxy 
laminates  are  used  as  reinforcement.  Consider  a  2D  analysis  in  which  Ej^  is  the  major 
modulus  of  the  laminate  and  Ey  represents  the  modulus  perpendicular  to  the 
laminate.  In  this  case  the  shear  modulus  is  taken  as  that  for  a  unidirectional  layer  of 
boron/ epoxy.  Using  typical  properties  for  the  boron/ epoxy  and  making  a  comparison 
with  the  ID  d.e.  on  the  basis  of  E2  =  E^  gives  the  restdt  shown  in  Table  6  labelled  as 
'orthotropic'.  Comparison  of  the  shear  stresses  obtained  by  the  F.E.  and  ID  d.e.  show 
that  better  agreement  is  obtained  with  the  F.E.  midplane  values  than  those  obtained 
using  equation  (13).  Overall,  the  comparison  with  midplane  values  in  Table  6  has 
shown  agreement  to  witiiin  6%. 
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4.3  3D  bonded  joints 

A  3D  F.E.  diennal  analysis,  using  20  noded  brick  elements,  has  been  carried  out  for  the 
structure  defined  in  Fig.  5,  which  is  fuUy  restrained  at  the  right  hand  end.  Also 
bending  in  the  z  direction  has  been  restrained.  In  this  case  the  structure  is  subject  to  a 
uniform  temperature  of  100°  C.  The  resulting  shear  stresses  have  been  evaluated  in 
the  middle  plane  of  the  adhesive  at  locations  A  and  B  as  shown  in  Fig.  5.,  and  have 
been  found  to  reach  a  maximum  value  at  a  distance  of  approximately  0.2mm  away 
from  the  end  of  the  joint. 

Results  shown  in  Table  7  correspond  to  an  isotropic  material  but  with  a  modulus  equal 
to  that  of  a  boron/epoxy  laminate.  This  is  simply  to  make  a  comparison  with  the  ID 
d.e.,  and  orthotropic  results  to  be  considered  later.  The  coefficient  of  thermal 
expansion  is  taken  to  be  4.1x10-^  /°C  as  used  previously.  The  results  are  shown  in 
Table  7  in  which  the  number  of  elements  in  the  width  of  the  patch  has  been  increased 
to  demonstrate  convergence.  At  point  A  a  stationary  value  has  almost  been  achieved, 
however  convergence  has  yet  to  occur  at  point  B .  Values  of  shear  stress  at  point  A 
(47.54  MPa,  Table  7  are  considerably  higher  than  the  corresponding  result  (32.19  MPa, 
Tables)  from  the  ID  equation.  The  shear  stress  values  at  point  B  (34.59  Mpa,  Table  7) 
are  significantly  closer  to  the  ID  equation  result.  These  results  indicate  that  a 
significant  Poisson's  effect  is  evident  in  the  3D  case.  It  is  therefore  corrsidered 
appropriate  to  model  the  composite  as  an  orthotropic  material.  A  3D  orthotropic 
analysis  will  now  be  considered. 

4.3.1  3D  orthotropic  joint 

The  material  properties  of  the  basic  unidirectional  laminate  of  boron/  epoxy  used  in 
the  3D  analysis  are  given  in  Table  8.  For  a  thermal  analysis  of  a  multi-ply  laminate  it  is 
necessary  to  compute  the  effective  coefficients  of  thermal  expansion  in  the  material 
symmetry  axes.  Usually  manufacturers'  data  will  provide  the  longitudinal  (aj  and 
transverse  {a 2)  coefficient  of  thermal  expansion  for  a  unidirectional  laminate.  In  this 
case  the  longitudinal  coefficient  of  thermal  expansion  is  a  measure  of  the  fibre 
property  while  the  transverse  coefficient  is  a  measure  of  the  resin  property.  The 
effective  coefficients  of  thermal  expansion  in  the  material  S5nnmetry  axes  are  given  by: 


where: 


a2^=-t(a,Cos^O+a2Sin^0)  ...(20) 

n  *=1''  ' 

a2.=~t{(XiSin^O  +  cc2Cos^O]  ...(21) 

n  k=-r  ’ 

k  is  the  k‘^  ply 

0  is  the  angle  of  the  k  ply  to  the  material  symmetry  axes 
n  is  the  total  number  of  pHes. 
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Note  that  the  derivation  for  equations  (20)  and  (21)  assumes  a  balanced  lay-up  and  is 
restricted  to  in-plane  loadings.  It  is ,  however,  possible  to  derive  an  expression  to  cover 
all  other  cases.  For  a  3D  laminate  the  out  of  plane  coefficient  of  thermal  expansion  may 
be  taken  as  being  equal  to  the  value  in  ihe  unidirectional  lay-up. 

In  Table  9  a  comparison  has  been  made  between  the  coefficients  of  thermal  expansion 
for  a  unidirectional  laminate  and  a  lay-up  consisting  of  [Oj  ,±45,03]^ .  The  overall  effect 
of  the  Poisson^s  ratio  and  the  lay-up  is  to  produce  an  effective  coefficient  in  the 
longitudinal  direction  (aj  higher  than  the  unidirectional  lay-up  and  a  transverse 
coefficient  ( than  the  unidirectional  lay-up. 

The  results  shown  in  Table  10  are  for  a  configuration  in  which  the  patch  is 
homogeneous  orthotropic  and  representative  of  a  [O2, ±45,03  lay-up  in  which  the 

principal  material  symmetry  axis  is  in  the  x  direction.  For  the  purpose  of  comparison 
with  the  ID  equation,  the  coefficient  of  thermal  expansion  in  the  x  direction  is  taken  as 
4.1xlO^/oC  ,  while  the  coefficient  in  the  y  direction  is  taken  as  zero.  The  number  of 
elements  used  in  the  width  of  the  patch  have  been  found  to  have  a  smaU  effect  on  the 
results.  The  maximtun  shear  stresses  ocour  in  the  middle  of  the  patch  (  point  A  ) 
while  the  t^x  shear  stresses  on  the  edge  (point  B)  are  slightly  lower.  The  maximum 
stresses  are  sUghtly  higher  than  those  predicted  by  the  ID  lap-joint  equation.  Although 
not  tabulated,  an  analysis  has  been  carried  out  for  the  same  problem  as  previously 
considered  but  with  coefficients  of  thermal  expansion  of  6.24xl0-6/oC  in  the  x  direction 
while  the  coefficient  of  thermal  expansion  in  the  y  direction  is  taken  to  be  zero.  These 
results  gave  a  9.7%  reduction  in  shear  stresses  compared  to  the  results  shown  in  Table 
10.  This  intermediate  step  was  carried  out  to  allow  a  useful  comparison  to  be 
considered  next. 

Results  shown  in  Table  11  are  for  the  same  problem  as  previously  but  correspond  to 
coefficients  of  thermal  expansion  of  6.24x10^/*^^  and  16.96x1 in  the  x  and  y 
directions  respectively.  These  results  are  20%  lower  than  the  previous  results  in  Table 
10  and  together  with  the  reduction  of  9.7%  from  the  intermediate  step  indicate  that  an 
interaction  is  occurring  between  the  Poisson' s  ratio  effect,  and  the  two  coefficients  of 
thermal  expansion.  Clearly  a  3D  analysis  is  significantly  different  to  a  ID  analysis 
when  both  longitudinal  and  transverse  coefficients  of  expansion  exist,  however  the  ID 
lap-joint  equation  gives  very  close  agreement.  A  comparison  of  all  of  these  results  is 
seen  more  clearly  in  Fig.  6. 

4.3.2  Peel  stresses 

Usually  the  end  of  a  bonded  patch  is  tapered  such  that  the  thickness  of  the  boron 
reinforcement  is  small  and  for  the  case  of  boron  would  be  approximately  0.127  mm. 
This  will  resTilt  in  significantly  lower  peel  and  shear  stresses.  While  peel  stresses  have 
not  been  validated  here,  closed  form  solutions  are  available,  and  are  usually  derived  as 
a  function  of  the  shear  stress,  ref (4). 
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5.  Residual  stress  in  a  bonded  repair 

5.1  Bonding  process 

The  bonding  of  a  high  strength  patch  to  a  plate  involves  heating  of  the  patch  and 
adhesive,  and  localised  heating  of  the  plate.  This  bonding  process  may  he  carried  out 
at  a  temperature  of  between  80  -  120°C  depending  on  the  type  of  adhesive,  and  may 
involve  a  duration  of  several  hours.  When  cured  the  patch  and  plate  is  allowed  to  cool 
to  room  temperature.  As  a  result  of  the  cooling  and  the  (Afferent  coefficients  of 
expansion  between  the  patch  and  tiie  plate,  residual  stresses  will  occur  in  the  plate, 
patch  and  adhesive.  If  the  thermal  coefficient  of  expansion  of  the  plate  is  greater  than 
that  of  the  patch,  as  in  the  case  of  aluminium  in  comparison  to  a  boron  patch,  then 
tensile  residual  stresses  will  develop  in  the  plate  while  compressive  stresses  will  occur 
in  the  patch.  If,  however,  the  plate  is  rigidly  restrained  close  to  the  patch  and  the 
effective  coefficient  of  expansion  of  the  plate  is  lower  than  the  patch,  the  residual 
compressive  stresses  will  exist  in  the  plate  and  tensile  stresses  in  the  patch. 

It  should  be  noted  that  for  analysis  purposes  the  shear  modulus  of  the  adhesive  must 
be  taken  as  some  average  value.  This  is  because  the  shear  modulus  of  the  adhesive  in 
the  cured  state  varies  significantly  with  temperature. 

The  procedure  outlined  here  will  enable  the  calculation  of  residual  plate  and  patch 
stresses,  and  adhesive  shear  and  peel  stresses.  However  the  validation  presented  here 
will  only  involve  the  plate  and  patch  radial  stresses.  The  circular  plate  analysed  in 
section  2  has  been  reinforced  in  the  central  region  with  a  circular  patch  of  radius 
r  =  Rj.  For  the  F.E.  analysis  two  cases  have  been  considered.  Firstly,  the  patch  is 
superimposed  over  the  plate  elements  and  shares  the  same  nodes.  The  patch  is 
considered  to  be  welded  to  the  plate.  An  analytic  approach  will  be  used  to  validate 
these  results.  Secondly  the  adhesive  has  been  included  in  the  analysis  as  3D  elements 
while  the  plate  and  patch  are  included  as  2D  elements.  In  this  case  it  is  necessary  to 
model  the  adhesive  as  a  3D  element  to  allow  the  transmission  of  shear.  Since  the  plate 
is  reinforced  on  one  side  only,  bending  effects  wiU  occur.  As  a  result  it  is  necessary  to 
restrain  the  bending  effects  so  that  a  comparison  with  an  analytical  approach  can  be 
made. 


5.2  Analytical  expression  for  residual  stresses  in  a  bonded  repair 

In  a  bonded  repair  both  the  patch  and  plate  are  heated  to  the  cure  temperature  T .  In 
the  case  of  a  plate  m  which  the  edges  are  restrained,  the  heating  gives  rise  to  an  initial 
stress  in  the  plate  before  curing  of  the  adhesive  and  cooling  takes  place.  In  the  case  of 
the  patch  no  restraint  exists  and  the  initial  stress  is  negligible.  To  compute  the  initial 
stress  in  the  plate  we  have  the  following  expressions  for  displacement  and  radial 
stress,  from  ref  (2): 
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U  = 


aJ\  +  v) 


\  Trdr  +  Cr 


...(22) 


„^_Z^(l^]Trdr+-^C 

0  (I-  v) 


...(23) 


from  the  boundary  condition  that  m  =  0  at  r  =  obtain  the  integration  constant 

C .  From  this  the  general  form  for  the  radial  displacements  and  stresses  is: 


^  j _S.]° Trdr] 

Ro  U  0  Ro  0  J 


fir  1  +  V  ,1 


Specifically  for  r  =  the  displacement  is  given  by: 


31 

Rl 


2\ 


V  no/ 


...(24) 


.(25) 


...(26) 


Specifically  for  r<Rj  we  have  a  constant  state  of  stress  given  by: 


rr,  ,  (1  +  ^) 


{T,-To)  (  R] 


’^^^IH.RjIR, 


-( 4-0 


...(27) 


Equation  (27)  gives  the  required  initial  stress.  The  only  reason  that  an  expression  for 
displacements  have  been  derived  is  a  means  of  checking  the  F.E.  results.  Both 
armlytical  and  F.E.  results  are  contained  in  Table  12  for  a  circular  plate  of  Ri  =  162  mm. 
The  material  properties  correspond  to  those  for  the  plate  given  in  Table  1.  A 
comparison  of  F.E.  and  analytic  results  shows  agreement  to  5  significant  figures  for 
both  displacements  and  stresses. 

The  second  part  to  the  solution  of  residual  stresses  involves  the  analysis  of  the  plate 
and  patch.  We  require  the  state  of  stress  corresponding  to  a  cooling  of  the  plate  and 
patch.  If  we  start  at  a  temperature  O^C  corresponding  to  zero  initial  stress,  then  cooling 
to  a  negative  temperature  equal  in  magrutude  to  the  cure  temperature  (-T^  will  give 
the  required  stress  components.  A  summation  with  the  initial  stresses  will  give  the 
residudl  stress  state. 

In  this  analysis,  the  adhesive  will  not  be  considered,  and  the  patch  is  assumed  to  be 
rigidly  connected  to  the  plate.  Radial  stresses  and  displacements  will  be  obtained  for 
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the  skin  and  patch.  Consider  a  plate  shown  in  Fig.7,  which  has  properties  E^,  a^,  for 
r>Rj  and  overall  properties  E^,  a^, for  r<Rj.  Again  the  temperature  boundary 
conditions  are  given  by  Tj  being  constant  for  r<Rj  and  T=Tq  at  r  =  RQ.  The 
displacements  are  given  by: 


for  r>Ri 


for  r<Rj 


w,  =  j  Trdr  +  Qr  +  -^ 

r  R,  r 


Ug  =  — - - J  Trdr  +  C,r 

r  0 


The  stresses  are  given  by: 


...(28) 


for  r>Rj 


a,E,  ] 


^Trdr- 


(1-v^) 


Q(l+v)- 


C3(l-v) 


for  r<Rj 


g-Q  =  -  J  Trdr  +  -  C/1  +  v; 

r  0  (\-v  ) 


...(31) 


The  solution  of  these  equations  must  satisfy  the  following  conditions: 

(a)  The  displacement  m,  and  is  equal  air  =  Rj,  hence  from  equations  (28)  and  (29): 


C^r  +  ^=  J  Trdr  +  C,Rj 


...(32) 


^  n  J  I  i  \  / 

r  Rf  0 

(b)  Equilibrium  must  be  maintained  across  the  boundary  at  r  =  Rj,  using  equations 
(30)  and  (31): 


c/i  +  v;- 

a-v^) 


C,a-v)l  a^Ej,  'r 


:  XL  i  t  XL  t 

°  °  °  |7Vdr+  ■  °  C 

Rj  r  0^' 


...(33) 


(c)  Also  we  have  the  further  botmdary  condition  that  at  r  =  R^: 


u,  =  0  = 


X^O 
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At  this  stage  we  have  enough  information  for  the  evaluation  of  the  constants  Q ,  Cj 
and  Cj. 

It  is  more  convenient  to  have  the  equations  in  the  form  that  represent  a  patch  over  the 
skin  for  r<R,  as  shown  in  Fig.  7.  As  before,  the  skin  has  the  properties  E„  a, 
while  the  patch  has  the  properties  E^,  necessary  to  derive  an  expression 

for  in  terms  of  these  quantities.  From  equilibrium  consideratioirs  we  have : 

Eoto^^o  =  W,  +  E^t^)a,  =  E,t,a,  +  E^t^a^  -(35) 


hence: 

Ja,+saJ  (36) 

(1  +  ^) 

where: 

s  =  ...(37) 

E\ti 

The  expression  for  the  stress  state  in  the  plate  just  outside  the  patch  is  given  by 
equation  (30)  for  r  =  Rj  : 


(\-v)  (l+v)Rj 


...(38) 


We  will  now  derive  the  expressions  for  the  stress  state  in  the  plate  beneath  the  patch 
and  in  the  patch.  From  equation  (29)  with  r<Rj  ,  under  a  uniform  temperature,  the 

displacement  is  given  by: 


u  = 


2 


...(39) 


Since  the  displacement  is  the  same  in  both  the  patch  and  skin  we  have; 


(l  +  v)a,Tj  ^  ^ 

2  ‘  Rj 

(\  +  v)aJj  ^  ^ 

2  '  Rj 


...(40) 

...(41) 


where  the  displacement  u  corresponds  to  the  location  r  =  Rj 
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The  radial  stresses  for  the  plate  and  patch  are  given  by: 


0 


J  Tjvdr  + 


^  _  "2-^2 

0^2  —  2 


r  0 


J  TjVdr  + 


EA 

...(42) 

a-v) 

E'y  C'y 

2  2 

...(43) 

(l-v) 

Using  equations  (40),  (41),  (42)  and  (43)  we  have  the  expressions  for  die  radial  stresses 
in  the  plate  beneath  the  patch  and  in  the  patch: 


cr,  = 


Er 


'  a-v) 


cr,  = 


(\-v) 


-ajj+  — 

R, 


-(X2TJ  H - 

Rj 


.(44) 


...(45) 


A  circular  plate  with  a  patch  has  been  considered  in  which  the  region  Ri  =  162  mm  has 
been  cooled  down  by  a  temperature  of  -lOQoC  and  the  outer  boundary  has  been  set  to 
QOC.  The  material  properties  correspond  to  those  in  Table  4.  Results  for  the  analytical 
solution  of  equations  (38),  (44)  and  (45)  are  shown  in  Table  13  together  with  the  F.E. 
results.  Displacements  agree  to  five  significant  figures.  Values  of  stress  in  the  patch 
and  beneath  the  patch,  agree  to  five  significant  figures.  However  the  agreement  of 
stresses  in  the  plate  just  outside  the  patch  varies  from  four  significant  figures  at  Ro  = 
500  mm  to  two  significant  figures  at  Ro  =  5000  mm.  This  is  entirely  due  to  the 
'stretched'  F.E.  mesh  at  the  larger  radius. 

To  obtain  the  residual  stress  in  the  plate  beneath  the  patch  it  is  necessary  to  sum 
equations  (27)  and  (44),  but  with  Tj  =-Tj  in  equation  (44).  Hence  the  final  expression 
for  the  residual  stress  beneath  the  patch  is: 


r  1 


To^ 


(Tj-Ta)  r  R] 


2\si{Rj  I R, 


M  ^-0 

0)^0 


+— ;...(46) 
(\-v)  ^  ‘  r/  ^  ’ 


Since  the  initial  stress  in  the  patch  is  zero,  then  the  residual  stress  in  the  patch  is  given 
by  equation  (45),  but  again  with  7}  =  -7)  hence: 

a  =  —^ — (a^Tj-^—)  -(47) 

(l-v)  ^  ‘  r/  ^  ^ 

and  the  final  expression  for  the  residual  stress  just  outside  the  patch  is  given  by  the 
summation  of  equations  (27)  and  (38),  hence: 
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L  ,  0  +  ^) 

2  1^^  (1-v) 


iJr-Ta)  (  R] 


2\n{RjlRoy  Ro 


(^-0 


c 


(\-v)  (l  +  v)Rt] 


...(48) 


These  equations  now  give  the  residual  stress  in  terms  of  the  aire  temperature  T.  If,  for 
example,  the  cure  temperature  is  100°C  and  residual  stresses  are  required  at  a 
temperature  of  -40PC,  then  this  would  correspond  to  a  temperature  of  T  =  +140oC. 


5.3  Comparison  of  F.E.  and  analytic  results 

The  numerical  solution  of  these  equations  has  been  achieved  and  the  foUowing 
quantities  have  been  evaluated  for  the  comparison  with  F.E.  results: 

(a)  displacement  at  r  =  Rj 

(b)  residual  stress  just  outside  the  patch  at  r  =  Rj 

(c)  residual  stress  in  the  skin  beneath  the  patch  ( <r  j ) 

(d)  residual  stress  in  the  patch  (<72) 

In  a  similar  way,  a  summation  of  the  F.E.  results  has  been  carried  out  to  compute  the 
same  quantities  above. 

As  previously  mentioned,  the  solution  for  the  uutial  stresses  are  contained  in  Tablel.2. 
Also,  the  solution  corresponding  to  a  temperature  change  of  -lOQoC  is  shown  in  Table 
13.  Finally  the  residual  stresses,  which  are  a  summation  of  these  two  cases,  are  shown 
in  Table  14. 

The  complete  results  in  Table  14  correspond  to  a  range  of  outer  plate  radii,  .  In  the 
case  of  the  displacements  agreement  between  the  analytic  and  F.E.  results  is  to  five 
significant  figures.  In  the  case  of  the  residual  stress  in  the  plate,  just  outside  the  patch, 
agreement  varies  from  four  significant  figures  at  the  lower  outer  radii  to  three 
significant  figures  at  the  larger  outer  radii.  Values  for  the  residual  stress  in  the  plate 
beneath  the  patch,  and  in  the  patch  agree  to  four  significant  figures. 

So  far,  the  adhesive  has  not  been  considered  in  the  analysis.  However  F.E.  results  have 
been  obtained  in  which  the  patch  and  plate  have  been  coupled  using  3D  adhesive 
elements.  To  make  a  useful  comparison  with  the  previous  analytical  work  the  bending 
of  the  plate  has  had  to  be  restrained.  The  material  properties  of  the  adhesive 
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correspond  to  the  adhesive  parameters  given  in  Table  5.  A  comparison  of  these  results 
is  shown  in  Table  15.  The  introduction  of  the  adhesive  has  resulted  in  an  error  of  only 
2%  in  direct  stresses,  showing  that  the  use  of  a  closed  form  solution  is  sufficiently 
accurate  for  patch  design. 


6.  Discussion  on  thermal  stress  analysis  using  PAFEC 

6.1  Thermal  stresses 

As  mentioned  previously,  thermal  stress  analysis  is  carried  out  in  two  steps.  Firstly  a 
thermal  analysis  is  carried  out  using  'thermal'  elements.  Some  of  these  elements  are 
shown  in  Table  16,  and  the  complete  list  is  contained  in  the  relevant  PAFEC  manuals, 
ref (3).  In  Table  17  an  example  control  module  is  shown  for  a  steady  state  temperature 
solution  in  which  the  computed  temperatures  are  to  be  written  to  file  HOTl.  For  the 
thermal  part  of  the  analysis  no  restraints  or  loads  are  applied  to  the  structure. 
Boundary  conditions  are  most  easily  specified  as  temperatures  on  nodes  of  the 
structtue,  however  they  can  also  be  specified  as  heat  flux  or  a  combination  of  both  but 
not  each  at  the  same  node. 

The  temperature  module  contains  the  temperature  boimdary  conditions  at  specified 
nodes.  It  could  contain  the  temperature  in  the  entire  structure  if  required,  but  this  is 
not  so  convenient.  A  heat  flux  modide  can  also  be  used  if  necessary,  but  generally  it  is 
easier  to  specify  temperatures.  Any  serious  problem  usually  involves  the  specification 
of  temperature  at  a  large  munber  of  nodes.  As  a  result  use  of  the  PUPPIES  program, 
ref(5),  is  the  only  feasible  approach  to  this  problem.  Use  of  this  program  allows 
temperature  to  be  set,  for  example,  on  the  circumference  of  a  plate,  or  alternatively 
within  a  radius  of  the  centre  of  the  plate  such  as  that  analysed  in  section  3. 

It  should  be  noted  that  thermal  elements  may  be  orthotropic  in  that  the  thermal 
coefficient  of  expansion  or  conductivity  may  be  different  in  different  directions.  The 
PAFEC  thermal  solution  is  carried  out  in  phases  1,  2, 4,  6  and  7  and  it  is  not  necessary 
to  go  to  phase  9. 

After  the  thermal  solution  has  been  completed  then  the  structural  analysis  is  nm.  The 
equivalent  control  module  for  the  structural  analysis  is  shown  in  Table  17  and  instructs 
PAFEC  to  read  in  the  temperature  file  HOTl  from  the  thermal  analysis.  Clearly,  in  the 
structural  analysis  restraints  must  be  specified,  although  applied  loads  are  optional. 
Another  module  required  in  the  structural  lun  is  shown  in  Table  18.  This  module  is 
included  exactly  as  shown,  and  signifies  that  a  steady  state  temperature  analysis  is 
involved.  The  structural  analysis  must  have  the  same  number  of  nodes  and  elements 
as  in  the  thermal  solution,  otherwise  temperatures  wfil  be  assigned  to  the  wrong 
nodes. 
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6.2  Residual  thermal  stresses 

For  F.E.  analysis  there  are  several  methods  for  calculating  residual  stresses.  An  outhne 
of  these  methods  is  given  below: 

6.2.1  Method  1 

Consider  an  example  in  which  a  boron/ epoxy  patch  is  bonded  to  an  aircraft  wing. 
Usually  an  F.E.  analysis  is  carried  out  when  a  repair  is  being  made  to  primary 
structure.  The  F.E.  model  involves  the  patch  and  underlying  skin  and  sometimes 
includes  spar  elements.  The  most  convenient  analysis  would  involve  the  structural 
detail  only,  in  which  an  effective  coefficient  of  thermal  coefficient  of  expansion  is  used 
for  the  plate  being  reinforced.  Obviously,  the  reinforcement  would  be  analysed  with 
the  actual  coefficient  of  thermal  expansion.  The  effective  coefficient  of  thermal 
expansion  can  be  obtained  using  equation  (10),  which  is  derived  from  a  circular  plate 
with  fixed  (and  free)  edge  restraints  for  the  case  Rq/Rj^O.  Studies  on  cirailar  plates 
indicate  that  this  is  an  upper  boxmd,  and  hence  is  conservative.  In  this  case  the  residual 
stresses  are  obtained  with  a  total  of  two  PAFEC  runs. 

6.2.2  Method  2 

If  the  structural  model  of  the  repair  is  sufficiently  large  that  it  is  connected  to  large 
surrounding  structural  elements  such  as  spars,  then  the  residual  stresses  may  be 
obtained  as  carried  out  in  section  5.  This  procedure  involves  the  calculation  of  the 
initial  stress  for  skin  area  corresponding  to  the  cure  temperature.  The  second 
component  involves  the  complete  structure  subjected  to  a  temperature  equal  to  the 
cure  temperature  but  opposite  in  sign.  Again  the  residual  stresses  are  obtained  as  a 
summation  of  the  two  components.  In  this  method  the  coefficients  of  thermal 
expansion  are  the  real  coefficients  for  each  material  and  a  total  of  four  PAFEC  runs 
would  be  required. 

6.2.3  Comment  on  methods  1  and  2 

If  this  analysis  is  to  be  used  to  obtain  residual  adhesive  shear  stresses,  then  the  correct 
choice  of  modulus  for  the  adhesive  is  important.  The  modulus  of  the  adhesive  varies 
significantly  with  temperature  and  properties  based  on  the  'average'  from  cure  to 
operating  temperature  should  be  used. 

6.2.4  Methods 

This  method  is  that  in  which  the  variable  material  option  of  PAFEC  is  used  to  simulate 
the  entire  bonding  process  and  would  involve  only  two  PAFEC  nins.  However  it  has 
not  been  possible  to  rtm  this  option  at  present.  Until  this  error  in  PAFEC  is  corrected, 
methods  1  and  2  are  the  only  available  options. 
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7.  Conclusions 


In  this  report  the  validation  of  the  thermal  stress  analysis  capability  of  PAFEC,  for 
steady  state  temperature  problems,  has  been  carried  out.  The  F.E.  results  obtained  for 
the  heated  circular  plate  are  identical  to  tiiose  given  by  the  one  dimensional  analytical 
solution,  since  the  solution  is  radially  dependent  only.  Furthermore  the  results 
obtained  for  the  residual  stress  are  identical  to  the  analytical  solution.  In  the  case  of 
bonded  joints  small  differences  occur,  as  a  resxilt  of  the  comparison  of  a  two  (and 
three)  dimensional  F.E.  analysis  with  a  one  dimensional  analytical  solution.  Use  of  the 
procedures  developed  in  this  report  will  enable  very  accurate  assessments  of  the  stress 
states  in  bonded  joints. 
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Table  1.  Material  properties. 


Thickness 

(mm) 

Young's  modulus 
(MPa) 

Poisson's  ratio 

Plate  (Aluminium) 

1.0 

71016. 

0.3 

Patch  (Boron) 

0.5 

156000. 

0.3 

Adhesive  (FM300) 

0.254 

2273. 

0.35 

Table  2.  Thermal  properties. 


Coefficient  of  thermal  expansion 

(/OC) 

Conductivity 

(I/ms»C) 

Plate  -  Aluminium  (tti) 

23.E-06 

160. 

Patch  -  Unidir.  Boron  (0t2) 

4.1E-06 

- 

Table  3.  Comparison  of  temperatures  predicted  by  F.E.  and  analytical  expression. 


r  (mm) 

226.5 

291.0 

936.1 

2118.7 

3107.8 

3967.9 

Analytic 

90.227 

82.921 

48.853 

25.036 

13.865 

6.741 

F.E. 

90.227 

82.920 

48.854 

25.036 

13.865 

6.741 

Table  4.  Comparison  of  effective  coefficients  of  thermal  expansion  for  heated  circular  plate. 
(Ri=162.mm) 


Ro(mm) 

Analytic  tteff 

F.E.  tteff 

%  error 

500. 

9.01E-06 

9.014E-06 

0.0 

750. 

10.30E-06 

10.30E-06 

0.0 

1000. 

10.95E-06 

10.95E-06 

0.0 

2000. 

11.99E-06 

11.99E-06 

0.0 

3000. 

12.40E-06 

12.40E-06 

0.0 

5000. 

12.77E-06 

12.77E-06 

0.0 

Table  5.  Material  properties  for  bonded  joint.  (ai  =  23xlO^/’C,  ai  -  4.1x10^  f>C) 


ti  (mm) 

t2  (mm) 

ta  (mm) 

Vi 

V2 

Va 

El  (MPa) 

Ga  (MPa) 

3.6 

1.778 

0.254 

0.3 

0.3 

0.35 

71016. 

842. 

Table  6.  Maximum  adhesive  shear  stress  in  a  bonded  joint , midplane  values  in  parenthesis. 
(ai  =  23x10^  pc,  ai  =  4.1x10-^  pC) 


E2/Ei(MPa) 

D.E.  (MPa) 

F.E.  (MPa) 

%  Difference  (based  on 
F.E.  being  correct) 

1.0 

27.43 

28.84 

5.1 

(25.77) 

(-6.1) 

1.5 

30.06 

33.48 

11.4 

(30.15) 

±3 _ 

2.0 

31.70 

36.46 

15.0 

(33.20) 

J±Z) _ 

2.2  (Orthotropic) 

32.190 

35.49 

9.3 

(31.09) 

(3.5) 
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Table  1.  3D  isotropic  F.E.  results  for  maximum  adhesive  shear  stress,  corresponding  to  a 
uniform  temperature  oflOO°C,  (a\  =  23x10^  f'C,  ai  =  4.1x10-^ /^C) 


Number  of 
elements  across 
width 

2 

4 

6 

8 

Shear  stress  at  A 

48.43 

47.69 

47.54 

47.56 

Shear  stress  at  B 

36.95 

35.03 

34.59 

34.36 

Table  8.  Properties  of  unidirectional  Boron/epoxy,  ref(6) 


Ell  (MPa) 

Ell  /  E22 

V12 

G  (MPa) 

207000 

10.89 

.21 

4800 

Table  9.  Thermal  expansion  coefficients  for  unidirectional  Boron/epoxy,  and  values  computed 
for  a  Boron/epoxy  laminate. 


Coefficient  of  expansion  (/°C) 

Unidirectional 

[02/  ±45,  Oils 

ttL  (longitudinal) 

4.1x10^ 

6.24x10^ 

aT(tr2msverse) 

19.1x10^ 

16.96x10^ 

Table  10.  3D  orthotropic  F.E.  results  for  maximum  shear  stress,  corresponding  for  a  uniform 
temperature  of  100  °C,  (ai  =  13x10^  /°C,  at  =  4.1x10^  f^C,  or  -  0/°C) 


Number  of 
elements  across 
width 

2 

4 

6 

8 

Shear  stress  at  A 

37.59 

37.90 

37.95 

37.91 

Shear  stress  at  B 

36.74 

36.05 

35.37 

34.72 

Table  11.  3D  orthotropic  F.E.  results  for  maximum  shear  stress  for  a  uniform  temperature  of 
100  oQ  (ai  =  23x10-6 />C,  ai  =  6.24x10-6 />€ ,  aj  =  16.96x10-6 


Number  of 
elements  through 
width 

2 

4 

6 

8 

Shear  stress  at  A 

30.19 

30.28 

30.30 

30.26 

Shear  stress  at  B 

30.19 

30.00 

29.82 

29.60 

Table  12.  2D  circular  plate:  F.E.  results  for  initial  stress  and  displacement  at  cure  temperature 
of  100  oQ  ai  =  23x10-6  joq  (analytic  values  in  parenthesis) 


Ro 

Displacement  at  r=Ri=162mm 

Initial  stress  at  r  <  Ri  =  162  mm 

(mm) 

(mm) 

(MPa) 

500 

0.14602 

-141.89 

(0.14602) 

(-141.89) 

750 

0.16686 

-128.84 

(0.16686) 

(-128.84) 

1000 

0.17741 

-122.24 

(0.17741) 

(-122.24) 

2000 

0.19433 

-111.64 

(0.19432) 

(-111.64) 

3000 

0.20082 

-107.57 

(0.20082) 

(-107.57) 

5000 

0.20692 

-103.76 

(0.20692) 

(-103.76) 
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Table  13.  2D  isotropic  circular  plate  with  patch.  Displacements  and  stresses  corresponding  to  a 
temperature  of  -100  oQ  ai  =  23x10^  pC,  ai  =  LlxlO-^  /C  (analytic  values  in 
parenthesis) 


Ro 

(mm) 

Displacement  at 
edge  of  patch 
(mm) 

Stress  just  outside 

patch 

(MPa) 

Stress  in  skin 
beneath  patch 
(MPa) 

Stress  in  patch 
(MPa) 

500 

-0.11499 

127.93 

161.33 

-66.813 

f-0.11499) 

(127.92) 

(161.33) 

(-66.813) 

750 

-0.12618 

113.22 

154.32 

-82.213 

f-0.12618) 

(113.21) 

(154.32) 

(-82.213) 

1000 

-0.13189 

105.72 

150.74 

-90.067 

f-0.13189) 

(105.71) 

(150.74) 

(-90.066) 

2000 

-0.14125 

93.476 

144.88 

-102.94 

(-0.14125) 

(93.409) 

(144.88) 

(-102.94) 

3000 

-0.14493 

88.621 

14Z57 

-108.01 

(-0.14493) 

(88.507) 

(142.56) 

(-108.01) 

5000 

-0.14843 

84.082 

140.38 

-112.82 

(-0.14843) 

(83.975) 

(140.38) 

(-112.82) 

Table  14.  2D  isotropic  circular  plate  with  patch.  Residual  displacements  and  stresses, 
ai  =  23x10^  />C,  ai  =  4.1x10^  /°C  (analytic  values  in  parenthesis) 


Ro 

(mm) 

Displacement  at 
edge  of  patch 
(mm) 

Residual  stress  just 
outside  patch 
(MPa) 

Residual  stress  in 
skin  beneath  patch 
(MPa) 

Residual  stress  in 

patch 

(MPa) 

500 

-0.11499 

-13.96 

19.44 

-66.813 

(-0.11499) 

(-13.97) 

(19.44) 

(-66.813) 

750 

-0.12618 

-15.62 

25.48 

-82.213 

(-0.12618) 

(-15.63) 

(25.48) 

(-82.213) 

1000 

-0.13189 

-16.52 

28.50 

-90.067 

(-0.13189) 

(-16.53) 

(28.50) 

(-90.066) 

2000 

-0.14125 

-18.164 

33.24 

-102.94 

(-0.14125) 

(-18.23) 

(33.24) 

(-102.94) 

3000 

-0.14493 

-18.949 

35.00 

-108.01 

(-0.14493) 

(-19.06) 

(34.99) 

(-108.01) 

5000 

-0.14843 

-19.678 

36.62 

-112.82 

(-0.14843) 

(-19.79) 

(36.62) 

(-112.82) 

Table  15.  2D  isotropic  circular  plate,  patch,  and  3D  adhesive  elements.  Residual  displacements 
and  stresses,  ai  =  23x10^  pC,  ai  =  4.1x10^  pC  (analytic  values  in  parenthesis) 


Ro 

Displacement  at 

Residual  stress  just 

Residual  stress  in 

Residual  stress  in 

(mm) 

edge  of  patch 

outside  patch 

skin  beneath  patch 

patch 

(mm) 

(MPa) 

(MPa) 

(MPa) 

500 

-0.11551 

(-0.11499) 

128.31 

(127.92) 

161.91 

(161.33) 

-65.532 

(-66.813) 

Table  16.  Corresponding  elements  pr  thermal  and  stress  analysis 


Element 

Thermal 

Structural 

2  noded  flange 

39800 

34400 

6  noded  triangle 

39110  (39115) 

36110  (36115) 

8  noded  quadrilateral 

39210  (39215) 

36210  (36215) 

20  noded  brick 

39710  (39715) 

37110  (37115) 
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Table  17.  Example  CONTROL  MODULES 


Thermal  Analysis 

Structural  Analysis 

CONTROL 

CALC.STEADY.TEMPS 

SAVE.TEMPS.TO.HOTl 

FULLCONTROL 

TOLERANCE=10E-3 

PHASE=1, 24,6,7 

CONTROL.END 

CONTROL 

READ.  TEMPS.FR0M.H0T1 

FULLCONTROL 

TOLERANCE=10E-3 

PHASE=1,2,4,6,7,9 

CONTROL.END 

Table  18.  Additional  module  required  for  structural  analysis 

TIMES.FOR.THERMAL.STRESS.CALCULATION 

TIME 

0  _ _ 
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Deformation  of  adhesive 

Figure  3.  Definition  of  parameters  for  ID  bonded  joint. 


Figure  4.  Definition  of  parameters  for  2D  bonded  joint. 
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Boron  patch 

T 

7  t 

_ d± 

L=1000mm,  Lp  =450mm,  W=450mm 
Figure  5.  Definition  of  parameters  for  3D  analysis  of  bonded  joint. 


— —  3D  F.E.  orthotropic  aj=0 


—  Lap-joint  d.e. 

—  3D  F.E  orthotropic 


Figure  6.  Comparison  of  adhesive  shear  stress  for  a  range  of  cases. 
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